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One of the most interesting aspects of graphene is the tied relation between structural and elec- 
tronic properties. The observation of ripples in the graphene samples both free standing and on 
a substrate has given rise to a very active investigation around the membrane-like properties of 
graphene and the origin of the ripples remains as one of the most interesting open problems in 
the system. The interplay of structural and electronic properties is successfully described by the 
modelling of curvature and elastic deformations by fictitious gauge fields that have become an ex- 
perimental reality after the suggestion that Landau levels can form associated to strain in graphene 
and the subsequent experimental confirmation. Here we propose a device to detect microstresses in 
graphene based on a scanning-tunneling-microscopy setup able to measure Aharonov-Bohm inter- 
ferences at the nanometer scale. The interferences to be observed in the local density of states are 
created by the fictitious magnetic field associated to elastic deformations of the sample. 

PACS numbers: 



I. LATTICE DEFORMATIONS AND FICTITIOUS MAGNETIC FIELDS IN GRAPHENE. 

Since graphite monolayers started to be isolated in a controlled way [U |2] , graphene has been an optimal play- 
ground to test the most exciting ideas in condensed matter A great deal of attention was initially paid to its 
striking electronic properties, but it was soon realised that the structural and mechanical properties can be even more 
interesting both from a fundamental point of view as well as for applications j4j ■ The modelling of curvature by gauge 
fields in graphene was suggested in the early publications associated to topological defects needed to form the fullerene 
structures [51 [3] . The main idea was that the phase acquired by an electron circling a pentagonal defect is the same as 
that arising when circling a solenoid with the appropriate magnetic flux in analogy with the Aharonov-Bohm effect. 
The fictitious magnetic fields were later applied to model the observed ripples [3-[9] and elastic deformations p!0Ul5j . 
The state of the art and an updated list of references can be found in [TB]. 

A turn of screw took place when the fictitious magnetic field became an experimental reality after the suggestion 
that Landau levels can form associated to strain in graphene |17l I18j and their subsequent observation in 19]. The 
formation of Landau levels requires very high values of the fictitious magnetic fields and hence very strong deformations 
of the samples. Fields of up to 300 Tesla were estimated for the observed nano bubbles in (TH]. In the present work, 
we consider the opposite limit where the fields are small and the electronic excitations can still be described in terms 
of plane waves (rather than Landau levels). As we show below, non-trivial amusing effects can also be generated 
by fictitious gauge fields in this limit. Specifically, we discuss the realisation of the Aharonov-Bohm (AB) effect 
via deformation fields and propose a simple scanning-tunneling-microscopy (STM) setup that can make use of this 
realisation to detect stresses. 



II. AHARONOV-BOHM INTERFERENCES FROM FICTITIOUS GAUGE FIELDS. 

The AB effect [20] is one of the hallmarks of quantum physics and has continued awaking interest along the years 
[SUES], partly due to the difficulty of finding appropriate experimental realizations. In our setup, the fictitious gauge 
fields arise from strained (or curved) portions of the graphene sheet and the AB interference manifests in the form of 
oscillations in the local density of states which can be detected by a STM. A simple way to generate this interference is 
by introducing scattering centres and thus engineering the scattering loop depicted in Fig. [T^. The electrons injected 
by the STM tip can follow this loop either clockwise or anticlockwise, which results in different phases due to the 
fictitious fields inside the loop. This setup does not require the two physical paths of conventional AB devices, and is 
analogous to the one recently proposed in [53] to detect nanoscale variations of magnetic fiuxes. In what follows we 
will analyse the different aspects of the fictitious-field case and show the experimental feasibility of our proposal. 

Time-reversal symmetry breaking is a key point in the magnetic AB effect. This gives rise to the additional AB 
phase in the electron wave function that, for closed trajectories, is proportional to the magnetic flux traversing the area 
enclosed by the loop. Thus, closed trajectories travelled clockwise and anticlockwise produce phases with opposite 
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FIG. 1: Basic sketch of the proposed device, a. The two scatterers (red balls) and the STM tip design a semiclassical closed 
path that electrons can follow clockwise or anticlockwise (red and blue lines in the figure). If the region enclosed by the path 
is strained or curved, the fictitious magnetic field associated to the deformation will give rise to a characteristic interference 
pattern in the LDOS arising from the asymmetry of the clockwise and anticlockwise paths induced by the fictitious field as 
described in the text (see also fig. 2). b Typical LDOS map to be observed in the STM experiment when the region enclosed 
is perfectly flat. 



sign, which can further be exploited as described in 23J, for example, to build a STM AB interferometer. In the 
case of clastically deformed or curved graphene, however, at first glance it seems that the fictitious fields will give no 
net AB effect because these deformations do not break time-reversal symmetry. The symmetry is preserved by the 
fictitious magnetic fields because these fields couple with opposite signs to the two different Dirac cones of graphene 
[6] (see the supporting information). In the following we show that, contrary to this initial expectation, a non-zero 
AB phase does arise from the elastic deformation fields in graphene. 

To our purposes it is convenient to write the electron wave function as '1'+ = {'^Kj'^'k') where 'i'K{K') is the 
bi-spinor corresponding the Fermi point K(K'). As discussed in the supplementary information, a deformation of the 
graphene lattice gives rise to a fictitious gauge field 

^ ^ / - Uyy\ 

where k ~ 3 nm^^, $0 is the flux quantum and Uij is the strain tensor which can be written as Uij = -^[djUi + 
diUj + {dih){djh)] in terms of the in-plane and out-of-plane displacements u and h respectively [IS]. The graphene 

wavefunction then accumulates a phase factor that can be written as e~"*o^^ for the electrons travelling along a 
closed path C, where T3 — diag(l, —1) is the third Pauli matrix and $ = A{\) ■ d\ is the flux of the flctitious field 
that passes through the corresponding area. Consider now the closed scattering path shown in Fig. [T^ and let a_ 
and a+ be the scattering amplitudes for electrons in the state ^0 that travel this path clockwise and anticlockwise 

respectively. These amplitudes can be expressed as a± = e "^*o^^ao, where ao is the scattering amplitude in the 
absence of field. The scattered wavefunction then will be such that (^E'ol^'s) = 2ao cos (7r$/"I>o). In plain words, 
the key issue is that the opposite gauge potential couples to different types of particles (electrons at each valley) and 
the net effect of the opposite trajectories adds up which is the same reason that makes it possible the formation of 
Landau levels discussed in [T71 [T5] . This physics is also very similar to the one discussed in [TOl HJ in relation with 
the absence of weak localisation in graphene. This mismatch between the electron wave functions will give rise to a 
distinctive signal in the local density of states (LDOS) that can be measured by an STM tip as follows. 

The LDOS can be expressed as N{t,(jj) = — f Im Tr G(r,r;w) in terms of the retarded Green's function G. In 
the presence of a (real or fictitious) gauge field A, it follows from the previous discussion that the electron Green's 
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FIG. 2: Strain-induced pseudo magnetic fields and AB interferences in tlie LDOS. a. Distribution of pseudo magnetic field 
generated by the in-plane strain of equation (4) discussed in the text. b. Interference correction to the LDOS at energy u = 0.4 
eV after subtraction of the A'^u=o contribution, as given by eq. (3). c and d represent the same for the case of the out-of-plane 
gaussian height profile of equation (5) discussed in the text. 



function for a homogeneous graphene sheet factorizes as 

Go(r - r', c.) - exp A(l) • dl J Goo(: 



r - r ,0; 



(2) 



in the semiclassical approximation, where 6*00(1" — r') is the Green's function for zero gauge field and the line integral 
goes along the straight line connecting r and r'. When the graphene sheet contains two scattering centres as depicted 
in Fig. [1^ the corresponding the LDOS turns out to be 



N{r,oj) = NA=o{r,uj) + Nioop{r,oj) {cos[TT<^{r)/<^o] ~ 1} . 



(3) 



where A'^^o represents the total LDOS in the absence of deformations, A^ioop the interference of all the scattering paths 
passing through the two impurities and enclosing a finite area in the absence of deformations, and $ the flux of the 
pseudo-magnetic field induced by the deformations through this area. Here the cos appears due to the AB interference 
associated with the deformation of the graphene sheet. There is a correspondence between this interference and the 
one described in |23| for the real magnetic case and non-relativistic electron systems. This correspondence, however, 
is rather non-trivial due to the matrix structure of the electron Green's function in graphene. The details of the 
calculation that gives rise to the expression ([3]) in graphene are given in the supplementary information. 



III. EXPERIMENTAL PROPOSALS. 



The fact that there is a real AB effect in the LDOS due to fictitious gauge fields that appear in deformed graphene 
is illustrated in Fig. [2j In a) we consider that the deformation of the graphene sheet corresponds to in-plane 
displacements (in polar coordinates): 

Ur = uor'^ sin39, ug = uor"^ cos 36 (4) 

as these described in [TH] but localised in a region of width a = 6 nm [uq = uoo/f ^ exp(— r^/2cr^) with uqq a constant]. 

The fictitious magnetic field configuration corresponding to this situation is shown in Fig. 2 a. Fig. 2 b. shows the 
interference pattern obtained at an energy of uj — 0.4 eV and Moo = 0.1 nm after subtraction of the Nu=o contribution 
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FIG. 3: Proposed STM interferometer as can be used to measure the strain in suspended graphene. In the figure the two 
impurities (red balls) are located in the sample at opposite extreme from the tip. The tip can stay at a fixed position while 
the fictitious magnetic field is changed by applying pressure either mechanically or with an electric field as described in [25j . 
The variation of the LDOS measured can be directly related to the mechanical deformation of the sample. The device can also 
be mounted on a graphene-sealed micro chamber with a small hole. When gas is pumped inside the chamber the graphene 
covering the hole will infiate as a balloon as described in i26j. The STM oscillations will be a direct measure of the tension 
generated by pressure in the graphene sample. 



for the given geometry. Figs. 2 c. and d. show the same computation for a pure out-of-plane deformation of the 
graphene sheet corresponding to a gaussian bump of profile [TH [M] 

/i(r) = ^exp(-rV2cr^), (5) 

We have done the calculation with cr = 6 nm, A = 0.5 nni. These two cases illustrate rather generically what can be 
expected. 

The physics discussed above opens the door to build STM devices to probe directly local stresses, and more generally 
mechanical properties, alternative to more indirect procedures as those using topographic measurements |27j or Raman 
scattering [28 . It can be of use in a variety of samples to characterise strain engineering designs coming from the 
substrates or growth conditions [29 . This can be illustrated with two point scatterers as before. A reference LDOS 
can be acquired by scanning with the STM tip and calibrate the device before producing the stresses in the material 
either thermally [30], mechanically [31], or from the corresponding substrate [32]. Then, by moving the STM tip 
along a line parallel to the scatterers, in such a way that the area of the scattering path is conserved, a map of the 
corresponding deformations can be obtained by subtracting the reference LDOS. Alternatively the setup can be used 
to characterise the samples by subtracting the Moop in Eq. ([s]) computed for an ideal clean, flat sample with the 
given positions of the scatterers. 

The suspended graphene flakes discussed in [3S] represent a system particularly suited for exploiting the AB effect 
due to fictitious gauge fields (see Fig. [s]). These samples are under tension and the fictitious magnetic field is maximal 
at the edges. This field has been estimated to be of a few Tesla in sub- micron wide flakes [33] . The STM tip 
and the scattering centres can be set in the supported (flat) region at opposite extremes from the suspended region 
whose mechanical deformation can be modulated by the electric force between the flake and the gate or by applying 
mechanical pressure. If the STM tip is let fixed at a given position, the cos $ oscillations of the LDOS signal gives a 
direct measure of the varying strain. 

Alternatively the STM setup can be mounted on the graphene-sealed micro chamber described in [26 . There, a 
bump in the graphene sheet can be created in a controlled way by pumping gas in the chamber, which will generate a 
fictitious magnetic field localised in the region inside the bump. Scattering centres then can be engineered in such a 
way that the fictitious field associated with the bump and the stresses generated here can be conveniently extracted 
from the LDOS. In particular, by placing the scatterers in the flat region away from the bump and performing STM 
scans parallel to them the area of the scattering loop is kept constant. Thus, a flrst scan done before creating the 
bump will show typical Friedel oscillations due to the impurities. Then, when a second scan is done along the same 
line after the gas in pumped and the bump forms, a different interference pattern will be observed in the dl/dV curves 
from which the distribution of pseudo-magnetic fltrx can be easily extracted. 

We note that, since the AB effect is due to the phase of the wave function, other candidates to produce AB 
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interferences in graphene are topological defects generated by substitution of a hexagonal carbon ring by a different 
polygon. These defects are very frequent in carbon nanotubes and other graphitic structures, pentagons and heptagons 
being energetically the most favourable. The fictitious fields created by these type of defects are non-abelian in the 
valley and sub-lattice spaces since they break the sub-lattice symmetry as described in In this context, non 
abelian means that the vector potential has an additional matrix structure acting in the valley space. The corrections 
to the density of states associated with these topological defects were analysed within the covariant formalism in [34] . 
The analysis of section 2 can be easily extended to deal with these fields. Taking into account that the square of extra 
matrix associated to the valleys (also a Pauli matrix) is the identity, our STM setup can also be used to probe the 
topological defects mentioned above. The phase interference will be proportional to the deficit angle of the defect and 
this possibility was already discussed in [35]. However, a possible drawback in the case of topological defects is that 
the effective magnetic field being strong and localised - they are in fact like solenoids - may act as an extra scattering 
centre spoiling the simplicity of the analysis performed in the present work. 

IV. PHYSICAL FEASIBILITY OF THE PROPOSAL. 

The physical feasibility of the proposed experiments relies on the capability to select and manipulate the appropriate 
scatterers and, most importantly, to ensure the validity of the semiclassical approximation. In the following we discuss 
this feasibility. 

Different sorts of scatterers can be considered for the setup. Vacancies, adsorbed atoms or even two passive STM 
tips can be used nowadays in a controlled way |36H38) . Inter-valley scattering will not suppress the effect but it should 
be minimised for the best performance of the proposed device. It will be suppressed if the range of the scattering 
potential is larger than the lattice constant but smaller than the electron wavelength. This range also allows to model 
the scattering potential as a i5 function in the continuum limit |39j . The best candidate for this would be adsorbed 
noble-gas atoms such as Ar or Xe, whose potential is intrinsically short-ranged. These atoms can be placed in the 
chosen positions for the experimental setup with the same STM tip (the process of single-atom manipulation with 
STM [40) is by now a well controlled procedure |41|). Since these impurities have lower migration barriers than 
covalently bond impurities, the experiment will be better controlled at low temperatures to ensure that the impurities 
will not move during the experiment [42j . Ionised impurities such as K with long-ranged Coulomb potentials can also 
be considered if a finite number of (gate-induced) carriers are added to screen their effects. 

The validity of the semiclassical approximation requires an electron wavelength much smaller than the magnetic 
length. The electron wavelength in graphene is A = 2itvf/E, where vp — 10^ m • s~^ and E is the energy measured 
from the Dirac point. We consider E = 0.4 — 2 eV ^ with the largest A ~ 10 nm for E = 0.4 eV . (These energies are 
routinely available in STM experiments [19]) . At higher energies in the case of graphene trigonal warping comes into 
play but this deviation of the band structure from circular shape is not a problem for the setup and can be easily 
taken care of. 

The distortions described in Eq. Q produce an effective magnetic field that can be safely estimated to be 5eff ~ 
k(woo/ct^)$o- This gives rise to an effective magnetic length ~ a / {kuqq)^^"^ . For (7 = 6 nm and Uqo 0.1 nm as 
considered before we have as > \E=OAeV^ which makes it possible to work within the semiclassical approximation as 
we did. Note also that these parameters make it possible to reach fluxes (j) ^ Bgjjcr^ such that (t>'r^/4)Q is of order 1. 
Similar numbers are obtained for the out-of-plane distortions described by Eq. ([5]) with A= 0.5 nm and cr=6 nm as 
we considered before. 



V. A SUMMARY 

As a summary, in this work we put together two very beautiful physical facts: the Aharonov-Bohm effect, and 
the physical reality of the fictitious magnetic fields associated to elastic deformations or curvature in graphene. The 
maturity of the graphene field and the experimental capability to produce and manipulate large samples of very high 
quality together with the technical advances in the STM devices makes it possible to perform a new and very original 
AB experiment. As a byproduct and based on the physics described, a device can be made to measure micro strains 
in graphene samples at nanometer distances using an STM interferometer. The cases provided in the paper are only 
illustrative examples of the possible devices and applications of the setup. 
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VII. SUPPORTING INFORMATION 

VIII. SUPPORTING INFORMATION 

A. Fictitious gauge fields and the continuum model for graphene. 

As it is well known by now, the 7r-electrons of graphene that are responsible for the electronic properties can be 
well described by a tight binding model in the Honeycomb lattice [43l |44] . The Honeycomb lattice can be seen as 
two interpenetrating triangular sub lattices. It has two points per unit cell. Due to this special topology, the Bloch 
function has two components representing the electronic probability at each sub lattice. The system at half filling 
instead of a Fermi line has six Fermi points located at the corners of the hexagonal Brillouin zone. Only two of them 
are inequivalent and can be chosen as Ki = (47r/3\/3, 0) and K2 ~ — Ki. Expanding the dispersion relation around 
one of the Fermi points, say Ki (fc^ — Kn + 5ki)^ one arrives to a continuum model for the low energy electronic 
excitations described by the equation 

where a ~ 1.42 A is the lattice constant and i ~ 3 eV is the tight binding From H/a in the limit a — > 0, t.a— const., 
one gets in the continuum limit the massless relativistic Dirac Hamiltonian 

Hi = hvpCT ■ k, (7) 

where Ui are the Pauli matrices and vp = ?>ta/2 c/300 is the Fermi velocity of the electrons. The same expansion 
around the other Fermi point gives rise to a time reversed Hamiltonian: H2 = hvp{—axkx + (Jyky). Notice that the 
two effective low energy Hamiltonians are related by the change cr 4-> ~cr* . 

The two Fermi points are referred to as "valleys" in a semiconductor language. In the absence of interactions or 
disorder mixing them, the valley quantum number is a degeneracy of the spectrum. We must note that time reversal 
is not an invariance of the effective hamiltonian around a given Fermi point but only of the full system with the two 
Fermi points. The scattering potentials and effective gauge fields invoked in this work may act on the valleys on a 
non-trivial manner. In the case of the scattering potentials, since the two Fermi points are far apart in momentum 
space (IK2 — Ki| ^ SX), to avoid intervalley scattering their range has to be larger than the lattice unit a, which is 
a reasonable assumption for the most common types of scatterers. 

Concerning the effective gauge fields associated either to elastic deformations or to intrinsic curvature of any origin, 
the crucial fact is the they couple to each Fermi point with opposite signs. The effective gauge fields can be obtained 
directly from the tight binding model if the derivation of the continuum model is done with arbitrary nearest neighbour 
hoppings ti different in the three directions [45] . The Hamiltonian around the Fermi point Ki is 

Hi = hvpo- ■ (-iV - A) , (8) 

where — ^ (^3 — 12) , Ay — ^ {t2 + — 2<i). The effective potential has opposite sign around the second Fermi 
point. The connection with elasticity theory is obtained by assuming that the anisotropy in the hoppings is due to 
atomic displacements of the positions u of the atoms: ti = t + 5ti, with dti ~ Ui. This defines the fictitious vector 
potential A given in Eq. 0. 

In the geometric approach to smooth curvature [12] the sign of the coupling of the effective vector fields depends 
on the definition of the spin connection and is determined by the product of gamma matrices in it. As the effective 
Hamiltonians around each of the Fermi points differ in the sign of only one of the gamma matrices the spin connection 
has opposite signs in the two valleys. In this approach the fictitious magnetic field is related directly to the intrinsic 
curvature of the sample. 

B. STM interferences in the density of states 

We now show how the quasiparticle interference pattern produced by two impurities in the LDOS is modified in 
the presence of a fictitious gauge field in a way analogous to the AB effect. The derivation mainly proceeds along the 
steps of Ref . [23 , though the fact that we have to deal with Dirac fcrmion Green's functions and the valley degree of 
freedom in graphene introduces additional subtleties that have to be worked out carefully. This is shown in detail in 
the following. 
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We first revise the case of zero strain, for one valley. The LDOS measured by the STM tip at the point r for the 
energy w from the Fermi level can be written as 



N{r,uj) 



ImtrG(r, r; w) 



(9) 



where G is the electron Green's function in graphene. For a clean flat graphene sample one gets the (zeroth order) 
Green's function 



Go(ri,r2; w) 



lUJ 



i/o(^|ri - r2|) + z^^^^-^i/i(c.|ri - ral 



(10) 



where Hq and Hi are Hankel functions [H] . The scatterers considered in the text can be associated with a scattering 
potential of the form U{r) = Uq [S{r — ri) + S{r — r2)]X, where I represents the identity in sub lattice and valley 
space. The Green's function can be expanded perturbatively in the scattering strength Uq as 



G(r,r) = Go(r,r) + J dr'Go(r, r')C/(r')Go(r', r) + 



(11) 



dropping the energy dependence since the scattering is elastic. This gives rise to the term Na=o{v,u!) in Eq. (|3|, 
whose precise form is unimportant for our purpose. 

Consider now the case of non-zero strain. Our purpose is to show that the contribution to the LDOS due to the 
gauge field can be recast in the simple interference-like term Moop (cos$ — 1) that appears in Eq. ([s]). As mentioned 
in the text, Moop is the part of Na=o generated by all diagrams associated with scattering paths in which the two 
impurities are involved (scattering loops). To see this, we compute the same perturbative series, but where the building 
blocks are now the Green's functions in the presence of the strain-induced gauge field. In the low-field regime, we can 
use a semiclassical approximation ^ for such a Green's function It is straightforward to show that terms involving 
one impurity only remain unchanged because these phase factors cancel. The loop terms are however modified. For 
example to second order we have 



(5G(2)(r, r) = U^tr [Go(r, ri)Go(ri, r2)Go(r2, r)e' W + Go(r, r2)Go(r2, ri)Go(ri, r)e~^ W 



(12) 



where obtained from the line integral J^^ + J^^ + , is just the flux of the pseudo magnetic fleld through the 
triangle. In the case of electrons in a normal metal, the Green's functions are scalars which moreover satisfy G(ri , r2) = 
G(r2,ri), so a cosine factor is obtained from this expression directly. The absence of these two properties makes the 
calculation in our case more involved. 

and then proceed to arbitrary order. We want to show 



We will first consider the second order case given by ( 12 



that the two terms in the expression ( 12 ) are actually equa 

G(ri, r2) = cr'^G(r2, ri)cr^ and the cyclic property of the trace to show that the second term can be rewritten as 



trGo (r , ri ) Gq (r2 , r) Gq (ri , r2 ) . 
We now commute the second and third Green's functions and obtain 

77$ 



(aside from the phase factors). For this, we first employ 

(13) 



(5G(2)(r,r) = U^tr 



Go(r,ri)Go(ri,r2)Go(r2,r)2cos ■ 



The commutator in this expression can be computed with (ITO 



Go(r,ri) [Go(r2,r),Go(ri,r2)] 



[Go(r2, r), Go(ri, r2)] ^ ^[a\ a^] ^'^^ '^^'['^^ '^^f Hi^r^ - r\)Hi{u\r^ ^ r^l). 
AvF |r2 - r||ri - r2| 

Since [a\ a^] — 2ie*%'^ for i,j = 1, 2, we have 

irGo(r, ri) [Go(r2, r), Go(ri, r2)] cx trGo(r, rOd^ ^ 0, 

so this proves that at second order the loop contribution is 

7r$ 



(5G(2)(r,r) = U^tr 



Go(r, ri)Go(ri,r2)Go(r2,r)2cos 



(14) 



(15) 



(16) 



(17) 



9 



Wc now consider a general term in the loop series, paired with the corresponding term where the path is traversed 
in the opposite direction. A generic term of this sort will contain multiple bouncing off a single impurity as well as 
back-and-forth scattering between the two impurities. The first type of terms are easily re summed in the T-matrix 
of the impurity (the energy index is also omitted) 



T = C/o + [/oG(rn,r„)f/o + --- = 



Uo 



l-C/oG(r„,r„)- 



(18) 



This expression is formally singular because of the real part of G'(rn, Tn). This has to be regularised, either by defining 
the real part in terms of the Kramers-Kronig transform of imaginary part, regularised with the bandwidth, or simply 
by an explicit computation in momentum space 



G(r„,r„) = J dPk 



— fc^ + ie 



- = w log 



to 



■ ITT 



(19) 



With A a momentum cutoff of the order of the lattice spacing. After this resummation, a generic term in the series 
will look like 

,„ ^, r -im .^^ .^^ 

^G(2™+2)(r,r) = irGo(r,ri)TiGo(ri,r2)T2[Go(r2,ri)TiGo(ri,r2)T2j Go(r2,r)e'^ + (1 o 2)6"'^. (20) 
Applying the same trick as before, and noting that the T-matrix is proportional to the identity, this turns into 
<5G('"'+'Hr,r) =irGo(r,ri)TiGo(ri,r2)r2[Go(r2,ri)TiGo(ri,r2)r2j Go(r2,r)e'*^ + 
trGoir, ri)TiGo(r2, r)Go(ri, r2)T2 [Go(r2, ri)T2Go(ri, r2)Ti 



e *o . 



So now we only need to commute the second Green's function to the right until it reaches the last place. As before, 
each step will produce a commutator proportional to a^. For example, the first will be 

trGoiv, rOTia^Ta [Go(r2, ri)T2Go(ri, r2)Ti 

Now the only non-zero traces in this expression contain and an even number of (T*(ra — ri,y with i = 1,2 and 
{fa — fb) the arguments of the Green's functions in the trace. All terms coming from the trace of such an expression 
will necessarily contain either vector products of proportional vectors (ri — r2)*(r2 — riY e'-^ , which vanish, or pairs of 
products of the type (r — ri)*(ri — r2)-'e'-' + (r — ri)'(r2 — riYe''^, which also vanish. This can be shown to happen 
for all steps in the commutation, so that the two paths are indeed equivalent except for the phase factor. 

Finally, a resummation of the back-and-forth terms can also be done, which leads to the final expression 



with 



and 



N{r,co) = Nu=oir,ij) + 7Vioop(r, w) (cos$ - 1) , 
iVioop(r, w) = Im trGo(r, ri)W^riGo(ri, r2)r2Go(r2, r) 



(21) 



(22) 



W 



l-TiGo(ri,r2)T2Go(r2,ri) 



(23) 



This calculation has been done for one valley. The contribution from the second one is obtained by changing 
and we get the same result because the cosine is even. 



